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ABSTRACT 



A new method for predicting the off design performance of a 
multistage compressor has been developed. This method replaces 
the usual routine numerical integration through the compressor by 
approximate analytical integration. The overall compressor perform- 
ance of the compressor is found as a perturbation of the design 
performance in terms of linear and quadratic expansions in two 
parameters, the deviation from the design flow rate and the deviation 
from the design speed. The procedure for this preliminary examina- 
tion has been simplified by assuming that all the single stage 
performance characteristics are identical and that the difference 
between static and stagnation conditions are negligible. 

The results of this method are quadratic equations, in terms of 
the two perturbation parameter s, for the temperature ratio, pressure 
ratio and adiabatic efficiency of the overall compressor. The 
amount of numerical work in finding the coefficients of these quadratic 
expansions is not negligible, however the method appears promising. 
Reasonably accurate results are expected for speeds varying as much 
as ten per cent from design and flow rates varying as much as four 
per cent from the design flow. 
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I. INTRODUCTION 

The usual technique for calculation of multistage compressor 
performance is numerical. Starting at the entrance with a prescribed 
flow rate and rotor speed the pressures, temperatures and flow 
coefficients are determined behind each stage in succession from the 
known performance characteristics of the individual stages. Although 
this method can be carried out rather rapidly on machines, it is 
desirable to have available a method that is more suitable for explor- 
atory investigations of the influence of stage characteristics. 

A new method is proposed which replaces the routine numerical 
integration through the compressor by approximate analytical integra- 
tion. The overall performance of the compressor can be found as a 
perturbation of the design performance in terms of linear and quadratic 
expansions in two parameters, the deviation from the design flow rate 
and the deviation from the design speed. The amount of numerical 
work in finding the coefficients in the expansion is not negligible, 
nevertheless the method appears promising. 

The investigation should be considered as preliminary. Cer- 
tain approximations, such as the assumptions that all single stage 
performance characteristics are identical and that the difference be- 
tween stagnation and static conditions are negligible, are not essential 
but simplify the procedure for the initial examination. 
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n. DEVELOPMENT OF EQUATIONS 



A. Work Coefficient 

The performance of multistage compressors is generally 
presented as plots of temperature and pressure ratios against dimen- 
sionless weight flow for given dimensionless compressor speed. The 
stage performance is most conveniently represented in the dimension- 
less form as a work coefficient (Y)# pressure coefficient (Y) ^.nd 
adiabatic efficiency {r\) against a flow coefficient (<[)). These 
dimensionless quantities are defined as: 




Ah^ is the increase in stagnation enthalpy per stage, c^ the specific 
heat of the gas at constant pressure, AT the increase in stagnation 
temperature per stage, Ap the increase in stagnation pressure per 
stage, and p the gas density. The compressor rotor reference 
speed is represented by U and the mean axial gas velocity by 
rfi is the mass flow of gas and A the cross sectional area of the 



flow annulus. 
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The multistage compressor temperature ratio, pressure ratio# 



and adiabatic efficiency will be determined from the single stage 
dimensionless characteristics. 

In general it can be said that the work coefficient Y is a 
function of the flow coefficient: 



The temperature rise (AT) in a single stage can thus be written as a 
function of the flow coefficient: 



A multistage compressor will be assumed to be made up of 
stages which all have the same characteristic. Such a compressor 



stages such that the temperature can be considered to be a continuous 
variable along the axis of the compressor. Let the continuous axial 
variable be replacing the sequence where n is the 

stage number from the entrance and N is the total number of stages 
of the compressor (0 ^ n i N). Using the variable ^ the tempera- 
ture rise can be written in the differential form: 



Y 





can be considered to be made up of an infinite number of infinitesimal 




The compressor then has an overall ’’length’’ of unity, and the 
function ^ (c|>) becomes: 



r(<t>)= N 



AJ V 
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The adiabatic compressor efficiency is normally defined as: 

= T. - T(°) 

T - rc”) 

where T is the actual temperature and is the isentropic final 

temperature for a given pressure rise. Introducing the isentropic 
relationship 




where and T/T(o) are the total pressure and temperature 

ratios^ p(o) and T(o) are the initial pressure and temperature and 
y is the ratio of specific heats of the gas. In a single stage the 
pressure and temperature increase can be represented by Ap and 
AT such that p = p(o) -f Ap and T = T(o) + AT. The stage 
efficiency thus becomes: 

(i + - I 



n = 



T 



Since Ap is small compared to p(o) the pressure term can be 
expanded to give: 



^ = 



ZiT 



^ 



Therefore in the limiting case, with Ap very small (as in the single 
stage): 
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^ ^P/p 



( 0 ) 



'VTi") 



or 



p(0) 



n 



H \ dT 

)^-i/ r(o) 



(1) 



This will also be the form of the limiting case as Ap — ► o so that the 
differential form can be written as: 



I (jp 

P di 



= ^ 



y ^ j_ JJL 

r-i) T ds 



(la) 



The adiabatic efficiency (r|) as defined here can be shown to 
be equivalent to the definition t] = limiting case. By 

the definitions of the work and pressure coefficients: 

Z) p 



V 

Y 



e(^Vz) 



Ap 



Cp AT 

UVz 



p Cp AT” 



But for Ap small equation (1) applies and: 

T' _ n(^)^ 



Y e 

Substituting the thermodynamic relations for a perfect gas: 

= R 



gives 



3L 

Y 
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B. Differential Equations for Flow Coefficient 

t has I 

t = 



The flow coefficient has been defined as: 

rh 



p U A 

Introducing the perfect gas relation, 

p = pRT 

rf,RT 



the flow coefficient becomes: 



P UA 



Taking the logarithm and differentiating with respect to ^ this 
becomes: 

J_ ii. _ JJL - 1 . L JA 



4 > 



or 



4 ' 



T d? 

T' 



p d€ A d^ 

A' 



r p A 

where the ( ) denotes differentiation with respect to the variable 



Substituting for ^ from equation (la) this reduces to: 



i: 

4 > 






~) 



y-i 



r 

T 



A' 

A 



( 2 ) 



This is an approximate relation since the differences between static 
and total temperature and pressure have been neglected in the 
development. 

At the compressor design point the flow coefficient, efficiency 
and temperature rise are assumed to be constants throughout and 
equal to the single stage parameters at the design point. Let the sub- 
script ( denote the design condition. Therefore at design 



conditions: 
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<j) = (j)^ = constant 



T| = r|^ = constant 



'ML 



( 3 ) 



I I ^ 

= y'= ^T= constant 

I o ^ O 



At the design point equation (3) can easily be integrated to give the 
temperature at any point in the compressor as: 



t: — 



u. 



z ^ -j- o. 



where is the inlet temperature. 

Dividing equation (3) by (4) gives: 



U." To 



Z 



To 



U." Z 



? - a 






(4) 



(5) 



where S is defined as: 
^o 






Uo' r. 



Equation (2) can now be solved to determine the area variation 
of the compressor: 



JL = 
A 



ri _i~i _n — —M. 
L i T (p. 



Letting the term r|( ^ly ) - 1 = f and noting that = 0 (since cj)^ 

is a constant) and substituting in equation (5): 

A' T 



A 



^ + ( 3 . 



( 6 ) 
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This area distribution is fixed by the design geometry. Substituting 
(6) into equation (2) the general differential equation becomes; 









_ f jr: + 

r 



(7) 



To eliminate the temperature (T) equation (7) is solved for T and 
substituted in the expression for T ; 



r = 



f T‘ 






_ Al 



r = 






r = 






•P r 



^ J 


fr 1 


— 


1 

t 'jji 


c/€ 


f. 

Lc+^. <p 




P. X 1 

4 + <P 






- 



Performing the differentiations and substituting the relations: 

<Jr dr d4> dr 



df 



d <|> d4> 

d f d 



d<p ^ 
df 4,' 



dS d(|> dC df 

there results the general differential equation in 

1J:\ f.(f-f.}r4> 









( 8 ) 



O 
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C. Solution for the Flow Coefficient 

A series expansion of each of the terms involved in the 
differential equation has been used to find a solution for the flow 
coefficient. From the single stage performance it is seen that both 
the efficiency function f and work coefficient T' are functions of 
the variable 

f = f (4>) r - T(<f>) 



These functions can thus be expressed as a Taylor series expansion 

about the design point. The stage efficiency at design is assumed to 

have its maximum value, hence the first derivative of the efficiency 

df 

function becomes zero ( - ^~ )^ = 0. The Taylor series expansions 
thus become: 






T(t)= T. 






(9) 



The flow coefficient will be represented by a series of the 



form: 



<j) - -h -I- e' + (10) 

The factor ^ represents a convenient small parameter to be 
replaced by unity in the final form. Substitution of the series (9) and 
(10) into equation (8) results in a series equation in powers of ^ . 
Equating terms in like powers of €. gives a series of linear 
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differential equations in the flow coefficients (n = 1 , Z , . , . n). 
Only the first two of these coefficients will be retained. The linear 
differential equations for these coefficients are: 






( 11 ) 




Note that the ( denotes differentiation with respect to the varible 

^ when used with the flow function ^ and denotes differentiation with 
respect to the variable <j) when used with f or ^ . The constants 

and (|)^ are all known from the single stage 

characteristic. 

To simplify the solution of equations (11) and (12) the following 
notations have been introduced: 









oC = 



r. 

-f ^ 






r. 

r. 



■ + 



f. 
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The coefficients of the left side of equation (12) become: 

i ^ f.) -f.H. ) — xf. 

Using this notation equations (11) and (12) become: 



/ " 

i>. + 



C< 



4 + ^. 



4>. “ o 



(Ha) 



4>"+ b4>V c Aiil +D 

* -c>r, T- (? +^.) 



(12a) 



Multiplying equation (11a) by the integrating factor (^ + p^) 

(4 •+"/3*) 4*. + +- /3*) -- O 

(?+/3.r 4>/ j - o 

Integrating: 

f? *• aT <i>.' = constant 

The constant can be found by introducing the boundary condition 
at ^ = 0; (j)j^ = <j)j^ (o) and therefore: 

Constant = (js^ (o) 

Therefore; 




• (lib) 



12 



Integrating between the limits 0 to 



- i (^) 



oC- I 

In terms of this becomes: 






(13) 



Equation 12a is solved in a similar manner. The integrating factor 
+ P^) is again applied to the left side of the equation and the 
solutions for and form (11b) and (13) substituted into the 

right side. On integrating once the equation is of the form: 

(})/-. E(?-(3.) ' 

The coefficients E , F, G, and H are a collection of all the constants 
in the equation not involved in the integration. On integrating the 
second term of the flow coefficient becomes: 



E -j-F ^ to' - '} 



( 14 ) 



The coefficients of equation (14) are defined as follows: 




# 



E 




c 

^ 0+1 



+ 



13 



iCX.-hz) 

] 0^ 

+ X, + I 




X* + 3 



-!> p.' 

G= 9- '♦'>' + ‘t>>> t,<»' + ^ 



( 15 ) 





(o) 



All terms conts-ining the variable ^ will be designated by lower case 
b's with subscripts: 





m 
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Introducing the coefficients (15) and (16) into equations (13) 
and (14) the first and second flow coefficient terms become: 



« (j)^(o) _ bj (j)/(o) (13a) 

(|)^ - (|)^(o) 4 - Q (|) (:o) -i-c (|)b) (|>(o)4- <J (|)^(o) 4- e (14a) 



where the coefficients of (14a) are defined as: 








(X.+I) 






VD ] 



d = D I ( by] 




Since a large number of coefficients will be used in the 
following equations, they will be represented by lower case letters 

with number subscripts. The letter will be used to indicate coeffi- 

2 ^ 2 
cients of a given term (a for 4 ) 2 ( 0 ); c for 4 ^j(o) 4 ^ 2 ( 0 ); d for 4^2 

and e for either 4^2 4^2^^^^* The number subscript refers to 

a given equation. The lower case b refers to terms containing the 

variable 



D. Flow Coefficient Initial Conditions 

Equations 13a and 14a can be used to determine the flow 
coefficient in terms of the initial conditions at ^ = 0; ( 4 > 2 (o)> 4 ^ 2 ^^^' 
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/ 



(j) 2 >(o) and 4)2 *(o)). The first stage of the compressor will be assumed 
to operate with the same characteristic as the single stage, therefore 
at ^ = 0^cj)(o) will be known. Let c|)(o) = c|)^ + €cl)j(o) with = 0 

for n > 1, thus <j>j(o) is prescribed and <t> 2 (o) = 0* For the initial 
condition at ^ = 0 equation (7) becomes: 




(17) 



By definition: 




T(^) = ® 




(18) 



The term 



Cp ® 



is similar to the term defined as , therefore let: 




( 19 ) 



Substituting from the expansions (9) and (10) gives: 





r(o) = 'i; -(-€?/ <{).(<’) 4- (p/of 

^(o) = f. + 4,/o)^ 




( 20 ) 
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Substituting (18) and the expansions (19) and (20) into equations (17) 
results in a series equation in powers of £ , Equating terms of like 
powers of £ gives the following expressions for the flow coefficient 
derivatives at the initial conditions: 













( 21 ) 


■(>>)- 


(%•/■ )' 





E. Flow Coefficient - Final Form 

Substituting the initial conditions (21) into the expressions for 
the flow coefficient terms (13a and 14a) results in a quadratic 
expression for the flow coefficient in terms of the velocity (or speed) 
perturbation pj, and the flow coefficient perturbation 4^|^(o) of 
the form: 

<p(4) + e -I- e"- 

The upper case letters will be used to represent coefficients of a 
given equation and the subscripts refer to the terms in the equation. 
The subscripts 1 and 2 will be used to represent coefficients of the 
terms <}>j(o) and respectively. 

Small deviations from the design point are accomplished by 
letting £ c|) 2 (o) and ^ be small. This can be accomplished by 
letting £ = 1 and using small perturbation of (t>j(o) and Pj^. 
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Letting £ = 1 the flow coefficient becomes; 



The coefficients of equation (22) are expressed in terms of the initial 
design conditions ) 9,) and the location along the length 

of the compressor These coefficients are most easily expressed 

in terms of the coefficients of all the previous equations involved. 

Express the initial conditions (21) as: 



<|) (o) = B, -1- §1 

cf>/o)= -t- 



(21a) 



where 




Substituting into equations (13a) and (14a) the coefficients of equation 
(22) become: 

A 2 ^ = — 

Alt =• a 

A,^ = £aB,Bi -t- c Bt — 

A, - ' - 

=» O! B| “1“ cBj + d — t>3 B|| 
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F. Temperature and Pressure Ratios 

The differential equation for the temperature rise through the 
compressor is integrated to determine the temperature ratios as 

z 

follows: JT U 



Tc?)— T(o) = 



U' 






Cp 



/ n^) 



Substituting the expansion of T' (4>) and dividing by T(o) = 0> the 
temperature ratio is: 

J . J. 



T(o) 



I + 






U T. 1 

Let £ = 1 and substitute — = rr the temperature ratio becomes: 

Cp 0 P 



r(§) 

Tio) 



I -I- + w.f'*- 



(23) 



The flow coefficient is defined as: 









U 



rh _ no R T~ 
()UA pUA 



The flow coefficient ratio at any station in the compressor becomes; 

4>(1) m R T(i) ^ p4») UM _ r(l) p(o A(o) 



4)(o) tfiRTC^) 'p(IMJtf) A^?) T(°) 1“^?) A(C) 

since rh and U are constant through the compressor. The pressure 
ratio is thus given by: 

TU) 4>(») A (o) 



> Co) 



T(o) A(€) 



( 24 ) 
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The temperature ratio is given by equation (23) and the flow coeffi- 
cients by equations (20) and (22). The area ratio is fixed by the 
geometry at the design condition given by equation (6) as: 



_L JA 
A 




Integrating gives: ^ e 

f M ^ ff 

•i A i 4^^- 



Ih 



ASH 

A(^'> 



AM 
A (5) 




The integrals in equation (2 3) are found by direct integration 
2 

of 4>j^, <j)j and <j >2 as given by equations (13) and (14). The results 
are again expressed in terms of the initial conditions and coefficients 
involving the design conditions and axial position (^) through the 
compressor. Thus: 

i / 

J <|>, d c = ^ <p,(0) -i- e, <P, 

0 

J t/'"’ + cl J. <!)/«> 

^ z z, 

fiJi- Oj <!><o)+ C3 +• dj -t- <£3 




m 
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In terms of the coefficients of equation (12) and the | variable 
coefficients (16) the coefficients of these equations are given by; 



e,- 









i.) ) (? + bj) ; d,- g 



r 

qJ-^^ fj 1 ,/t. . C-2P _ , P ^ 





CUX.i-i) 



2I> ~] 









^ ^ 2 . 
X>-hl 



-f- -h|3.) tj-J 




e 
^ ( 



Substituting these integrals into equation (23) and introducing 
the initial conditions as given by equations (21a) the temperature 
ratio also reduces to a quadratic expression in the velocity perturba- 
tion and the initial flow coefficient perturbation <j>j^(o). This 

equation is of the form: 
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-Ell 

T(») 



I 



-hj)^ + 4^. 4>/o)+D, t,<0> -Hl»»)j- 



(23a) 



The coefficients of equation (23a) are given by: 



!>.- 


To Ti 

r. 


e - ^ 






To' f 
To 


A (“iB. )+ 




D,- 


T.' 


(B,e,+ ?) 




To 




D,.- 


4( 


a,B,B, H- B, c, +- B. C,)- 


D. 

(3' 


D„- 


To' 

To 


(b; a, + 3, c. + d , 4-B„e,) + ^ (aX+ B, C,+ « 


0 



The maximum value of the temperature ratio at a given speed 
(Pj^ = constant) is found by setting the derivative of the temperature 
ratio with respect to <j)j^(o) equal to zero: 











The value of the flow coefficient for the maximum temperature ratio 
is thus: 




D.Z (3, 4- D, 

ZV„ 



The second derivative corresponds to the curvature of the temperature 



ratio curve and is: 
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r (jncU z-Po 
j <^.4 L J (3. 



and is independent of The temperature ratio curves can thus be 

completely determined by the coefficients of equation (23a). 

A similar parabolic curve can be found for the pressure ratio 
by substituting equations (20), (22), (2 3a) and (2 5) into equation (24), 
and dropping all the terms of and 4>j^(o) higher than the second 

order. The same procedure can be used to express the adiabatic 
efficiency r| by substituting the equations of pressure and tempera- 
ture ratio in the equation: 






u) \ » 



pCt) 



T(?) 

Tfo) 



- I 



The compressor efficiency (r|^) is determined by setting ^ = 1, 

The second order equations for the pressure ratio and efficiency thus 
reduce to the following form; 



>2 = ^» * -4- Q (3, C, (p(o)-t C„ (26) 



The coefficients of these equations are most easily expressed in terms 
of the coefficients of the previous equations and are: 
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Di 




Ai. 


. 




l-A. 












) 






-f- 


cr ™ 


IM- 


A 


Zl 


A.D 


^ -L 










<Po 


4>.C?^p.) 


l + p. 




r~ ^ 


2A,A, 


* At 


All 


1 


* Td 




- A, A - /AxAl 




<}>; 




4>. 






I2L * 


J 


E,r 


Br)' 


— — 


A, 

4>.‘ 


/All 






• Di A 1 

<p. J 




K 
K- 1 




- 


I>Z 


) c, 


K 

“ K-J 






K 

«-l 


x-i; 




m 

1 


EA- 


e.eB 
Y J 


+ 2.AA- 2),;*. 




K 

K-l 






£A- 


j 


+ dI - 


-A. 




K 

K-l 


x-i 

y 


{E„ 




1 

m 


} + Df 


-A 








Id 

k. V 


o^f.) 








K 


-(- 




) 




; ^ = 


e 

CQ„ 

II 


-0 



The coefficient Cj can be assumed to be zero since the 
maximum efficiency at design speed is known to occur at the single 
stage design flow rate. Setting the derivative of the efficiency with 
respect to the initial flow coefficient equal to zero: 

C. +2C„4>4=0 

the maximum efficiency is found to occur at an initial flow coefficient 
perturbation of: 
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cj)(o)=. — I _C. Ciz (3, 1 
L 2 Cn J 



Since the maximum adiabatic efficiency at design speed = 0) must 

occur very close to the design flow, should be very small. 

The speed function for maximum efficiency is found by 

substituting the value of the flow coefficient perturbation ((()^(o)) for 
maximum efficiency into equation (2.^) which becomes: 

n - +(c.- p. -f^) 

Setting the derivative with respect to equal to zero the maximum 

efficiency is found to occur at a speed function of: 
n — 2 Ch 
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ni. APPLICATION 

The performance of a ten stage compressor has been computed 
to illustrate the use of the method. The single stage characteristics 
for the research compressor as described in reference 1 have been 
used in this analysis. The single stage characteristics for this 
compressor are shown in Figure 1 in terms of the parameters ^ 
and f. The values of the design point constants for this example are; 

r = 1.95, r' = -7.25, r'" =- 37 . 5 , f = 2 . 01 , f'^ = - ns 

o o o ' o o 

and (|)^ = 0,44, The reference rotor speed for this example is 

= 1000 fps. The gas is air at standard temperature and pressure 

with S> = 520°R, C = . 24 Btu/lb°R and ^ = 1.4. The value 
o p 

of 3 and 2^ then become; 6 = 1. 603 and 2^ = - 3.288. The 

'^o o ^o o 

broken line curves in Figure 1 are the approximate curves represent- 
ing the quadratic expansions found by dropping all terms higher than 
the second order in equations (9). This approximation is seen to be 
reasonably accurate for flow coefficients from 0. 40 to 0. 50. 

The overall compressor performance is found by letting 
N = 10 and | = 1 . 0 in equations (23a), (25) and (26). The speed 

variable is computed from equation (19) which can be simplified to 
give: 




The resulting temperature ratio curves for the ten stage compressor 
are shown in Figure 2 as a function of the inlet flow coefficient. The 
pressure ratio and efficiency curves are shown in Figure 3. The 
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maximum efficiency is found to occur at Pj = 0. 363 which corre- 
sponds to U/U^ = 0.903, 

The operating point of each stage is found by varying ^ in ten 

equal increments and solving equation (22) for the flow coefficient. 

The operating points at design speed for the first, second, fifth and 

last stage of the ten stage compressor are shown in Figure 4. The 

stage flow coefficient (4>) is plotted as a function of the compressor 

inlet flow coefficient (cj) (o)). The stage efficiency (r|), work coeffi- 

/ 

cient (Y) pressure coefficient (1^) shown as functions of 

the stage flow coefficient (c|)). These curves indicate the range in 
which the quadratic approximation is reasonably accurate. The quad- 
ratic expansion of and f are shown in Figure 1 to be reasonably 

accurate for stage flow coefficients between 0. 40 and 0. 50. The tenth 
stage is seen to be outside this range at design speed for compressor 
inlet flow coefficients greater than 0. 455 and less than 0. 400. 

The last stage flow coefficient will have the greatest variation 
from the compressor inlet flow coefficient for off design speed 
conditions. The operating points for the tenth stage at varying speeds 
are shown in Figure 5. 

The operating range of the tenth stage indicates that the method 
should be reasonably accurate for a fairly wide range of compressor 
inlet flow coefficients at design speed. For varying speeds reasonable 
accuracy can be expected at speeds varying up to ten per cent off 
design with a similar flow coefficient deviation. 
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CONCLUDING REMARKS 

It has been shown that the off design characteristics of a multi- 
stage compressor can be predicted by a method using approximate 
analytical integration. The method can be expected to give reasonably 
accurate results at design speed for compressor inlet flow coefficients 
varying as much as four per cent from the design condition. At off 
design speeds the method should be accurate for speeds varying up to 
ten per cent from design with similar flow coefficient variation for each 
speed. A comparison with actual test results has not been made since 
none could be found for a compressor with all stages having identical 
characteristics. 

This method should be extended to include varying stage 
characteristics. This could be done by assuming a given single stage 
characteristic over a section of the compressor and different single 
stage characteristics over other sections. The approximate integra- 
tions could then be carried out over the various sections using the 
appropriate stage constants for that section. 
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Fig. 2 Ten Stage Compressor - Overall Temperature Ratio. 
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Fig* 3 Ten Stage Compressor - Overall Pressure Ratio 
and Adiabatic Efficiency. 
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